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Abstract. In the paper of Uwano [Czech. J. of Phys., vol.56, pp. 1311-1316 (2006)], 
a gradient system is found on the space of density matrices endowed with the quantum 
SLD Fisher metric (to be referred to as the quantum information space) that realizes 
a generalization of a gradient system on the space of multinomial distributions studied 
by Nakamura [Japan J. Indust. Appl. Math., vol.10, pp. 179-189 (1993)]. On motived 
by those papers, the present paper aims to construct a gradient system on the 
quantum information space that realizes the Karmarkar flow, the continuous limit 
of the Karmarkar projective scaling algorithm for linear programming. 
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1. Introduction 

There exist various excellent algorithms developed in engineering and systems science; 
as one of those, the Karmarkar projective scaling algorithm for linear programming is 
very famous, which solves linear programming problems in polynomial time (Karmarkar 
1984). 

In the middle of 1990's, Nakamura revealed an integrability behind several 
algorithms and governing equations arising in engineering and systems science; the 
Karmarkar flow emerging as a continuous version of the Karmarkar projective scaling 
algorithm (Karmarkar 1990, Nakamura 1994a), an algorithm solving the eigenvalue 
problem of anti-Hermitean matrices (Nakamura 1992), an averaged learning equation 
of the Hebb-type (Nakamura 1994b) and a gradient system relevant to the information 
geometry of multinomial distributions (Nakamura 1993). 

More than ten years after Nakamura's work (Nakamura 1993) on the gradient 
system on the space of multinomial distributions (GS-MD), one of the authors, 
Y. Uwano, encountered a very natural generalization of the GS-MD through a study 
on geometry and dynamics of a search algorithm for an ordered-tuple of multi-qubit 
states (Uwano 2006, Uwano et al 2007): The space of density matrices endowed with 
the quantum SLD Fisher metric which is referred to as the quantum information space 
(QIS) was constructed through a geometric reduction of the space of ordered-tuples of 
multi-qubit states. On the QIS, the gradient system associated with the negative von 
Neumann entropy (GS-NVNE) was studied: The matrix-form solution of the GS-NVNE 
was obtained explicitly, whose diagonal part was shown to describe the solution of the 
GS-MD. The GS-NVNE on the QIS is therefore understood as the natural generalization 
of the GS-MD. 

This encounter encourages the authors to seek other gradient systems on the QIS 
that realize certain dynamical systems or algorithms in engineering and systems science. 
The aim of the present paper is to construct a gradient system on the QIS that realizes 
the Karmarkar flow for linear programming. In what follows, the organization of this 
paper is outlined. 

In section 2, the general framework of gradient systems on the QIS is derived 
together with a review of the QIS: The construction of the QIS is accomplished by 
applying a geometric reduction method to the space of ordered tuples of multi-qubit 
states for a quantum search (Uwano 2006, Uwano et al 2007). After the review, the 
equation of motion is derived for arbitrary gradient systems on the QIS. Section 3 is 
the core part of the present paper, where the gradient system on the QIS realizing the 
Karmarkar flow (GS-QIS-K) is given explicitly. The Karmarkar flow mentioned here is 
the family of trajectories arising from the Karmarkar projective scaling algorithm for 
the unconstrained case (Karmarkar 1990, Nakamura 1993). A key to find the GS-QIS-K 
is to observe that the Riemannian structure of the canonical simplex, the underlying 
manifold, for the Karmarkar flow is in isometry to the Riemannian submanifold of the 
QIS consisting of diagonal matrices. Section 4 is devoted to the concluding remarks. 
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Throughout the present paper, differential geometric calculus works very effectively, for 
whose detail Appendices A and B are prepared. 

2. Gradient systems on the QIS 

2.1. Geometric setting-up for the QIS: Review 

Following Uwano (2006) and Uwano et al (2007), we review the geometric reduction 
method to construct the quantum information space (QIS) from the space of normalized 
orederd-tuples of multi-qubit states (STMQ). The natural Riemannian metric of the 
STMQ is shown to be reduced to the quantum SLD-Fisher metric of QIS. Those who 
are not so familiar to differential geometry may skip subsection 2.1, which is however 
deserves for tracing basic ideas for a series of works on gradient systems on the QIS by 
the author(s) (see also the closing remark in subsection 2.1). 

2.1.1. Reduction of the STQM to the QIS By M(2 n , to), we denote the Hilbert space 
of 2" x m complex matrices endowed with a natural Hermitean inner product 

($, $') = — tr^^') ($, $' G M(2 n , to)), (1) 
m 

where the superscript t will indicate the Hermitean conjugate operation from now on. 
The M(2 n ,m) is thought to describe the Hilbert space of ordered tuples of multi-qubit 
states, if every column, (f>k>) g C 2 ™ (j — 1, 2, . . . , to), of a matrix 

$ = (^,.-J (m) )eM(2 n ,m) (2) 

is understood to express a vector in the standard complex Hilbert space, C 2 ™, of n-qubit 
states (see Nielsen and Chuang 2000 for basic setting-up of quantum computation, e.g., 
multi-qubit states . . . etc). 

Remark 1 The set of q x r complex matrices is denoted by M(q,r) through this paper. 

We denote by M!(2 n ,m) the subset of M(2",m) consisting of $'s in M(2 n ,m) with the 
unit norm. Namely, 

Mi(2 n , m) = {$ G M(2 n , m) | ($, $) = 1}. (3) 

The space M 1 (2 n ,m) is understood as that of normalized ordered-tuples of multi-qubit 
states, which will be abbreviated to the STMQ. 

Let us consider the space ofmxm density matrices (see Nielsen and Chuang 2000), 

P m = {p £ M(to, to) j — p ; tr p — 1, p : positive semi definite}, (4) 

as the quotient space, say Mi(2™, m)/U(2 n ), of Mi(2",m) with respect to a natural left 
U(2 n ) action 

$ G M!(2 n , to) i — >■ (?$ G M!(2 n , to) ($ G M 1 (2", to), # G U(2™)), (5) 
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where U(2 n ) denotes the group of unitary matrices of degree 2 n . Indeed, if we define 
the map of M!(2 n , m) to P rn to be 

7r m : $ G Mi(2 ft ,m) .— > G P m , (6) 

m 

we see that n m is surjective and that vr m ($) = 7r m ($') holds true for $, $' G Mi(2 n ,m) 
if and only if there exists g G U(2 n ) subject to $ = g<&' . This shows 7r m (Mi(2 n , m)) = 
P m = Mi(2 n ,m)/U(2 n ). It follows from (TJD and © that the U(2 n ) action leaves relative- 
configuration among (fi^s (multi-qubit states) in an ordered tuple $ invariant, so that 
P m realizing Mi(2 n , m)/U(2 n ) can be referred to as the space of relative- configurations 
of multi-qubit states in ordered tuples (Uwano 2006 and Uwano et al 2007). 

2.1.2. The quantum SLD Fisher metric on the QIS To proceed differential calculus 
including differential equations, we have to consider a regular part of P m , which is 
realized as the mxm regular density matrices denoted by P m (Uwano 2006 and Uwano 
et al 2007), which can be referred to as the space of regular relative- configurations of 
multi-qubit states in ordered tuples. 

As the space oimxm regular density matrices, P m admits the quantum information 
space structure endowed with the quantum SLD Fisher metric denoted by ((■, -))® F . The 
quantum SLD Fisher metric is introduced as follows. 

Let us consider the space ofmxm traceless Hermitean matrices, 

T p P m = {S G M(m, m) | H f = 2, tr 2 = 0} , (7) 

as the tangent space of P m at p. The symmetric logarithmic derivative (SLD) at p G P m 
for H G T p P m is defined to provide the matrix C P (B.) G M(m,m) subject to 

i { P C P (E) + C p (E) p} = ~ (Ee T p P m ) . (8) 

The quantum SLD Fisher metric, denoted by ((•, -))® F ', is then defined to be 

((2, EX F = \tr [p(L p (E)L p (E>) + L P (E')L P (E))} (S, E' G T p P m ) , (9) 

(see Amari and Nagaoka 2000, Uwano et al 2007). 

We wish to express ((•, -)) QF explicitly. Let p G P m be expressed as 

p = hQh\ h G U(m) 

(10) 

9 = diag(0i, . . . , d m ) with tr6 = l, 6 k > (k = 1, 2, • - -m), 
where U(m) denotes the group ofmxm unitary matrices. Expressing E G T p P m as 

E = h X h ] (11) 
with h G U(m) in fflOl) . we obtain an explicit expression, 

(tf£ p {E)h) jk =—?—-Xik ti,k = l,2,..-m), (12) 

Uj + Uk 

of the SLD to E G T p P m . Putting pD]) - (IT2|l into ©, we have 



j,k=i °i + k 



'p 
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where E' G T p P m is expressed as 

E' = hx'hl (14) 

Surprisingly, the quantum SLD Fisher metric ((•, -)) QF thus defined turns out to be 
identical with the Riemannian metric, denoted by ((•, -)) R , that is determined along with 
the dimensional reduction of ir^-(P m ) (c M!(2 n ,m)) to P m through 7i m given by (Q. 
Indeed, under ( ITUl) . ( TTTj) and (fl4"l) . we have 

((S, 30)? = \ t f^f = tP. 3'))f (15) 

j,k=l 3 k 

(see Appendix A for the construction of ((•, -]) R and Uwano et al 2007). 

Theorem 2.1 (Uwano 2006, Uwano et al 2007) The quantum SLD Fisher metric 
((') ~))® F coincides with ((•, -)) R up to the constant multiple; 

((-, ~x F = 4 ((3 E X (3 s ' e w- ( 16 ) 

Throughout this paper, we will refer to the space of m x m regular density matrices, 
P m , endowed with the quantum SLD Fisher metric ((•, -))® F as the quantum information 
space (QIS) for brevity, that will be indicated also as the pair (P m , ((■, -))® F )- 

2.2. General framework of gradient systems on the QIS 

Now that the QIS is constructed as the Riemannian manifold with the metric ((•, ■))® F , 
we move on to study gradient systems on the QIS. The motive already mentioned in 
Sec. 1 for studying gradient systems on the QIS comes from the big similarity discovered 
by Uwano (2006) between the gradient system on the QIS with the potential equal to 
the negative von Neumann entropy and that on the information space of multinomial 
distributions studied by Nakamura (1993). 

Let F be a smooth real- valued function on P m . With the quantum SLD Fisher 
metric ((•, -))® F , the gradient vector, denoted by (gradF)(p), of F at p is defined to be 
the tangent vector of P m at p subject to 



(((gradF)(p),S')) p QF = (dF) p (S') = | 



F( 7 (t)) ( V 3GT p P m ), (17) 



t=o 



(see Kobayashi and Nomizu 1969). The ^(t) with a < t < b (a < < b) in (JTTJ) is a 
curve in P m subject to 

t G [a, b] ^ 7 (t) G p m , 7 (0) =p, -J =3. (18) 

Through a straight but quite long calculation, an explicit expression of the gradient 
vector (grad F) (p) is obtained in what follows. 

Let us introduce the real- valued variables, Xjk, yjk (1 < j < k < m) and zg 
{I — 1, 2, ■ ■ • , m), to express the entries of p G P m as 

Pjk = p k j= x jk + Wjk (l<J<k<m), pee = z e (£ = 1, 2, • • • , ra), (19) 
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where zfi are subject to 

m 

z e >0 (£=l,2,---,m) and JT^ = 1. (20) 

£=l 

For the complex variables p ab (1 < a < b < m), we introduce the partial differentiations, 

d 1/9 d \ d 

— i- 



dp, 



ab 



2 V&r 



ab 



dy a 



d d 



(1 < a < b < m) 



(21) 



dp ab 2 K -dx ab ' ' dy ab - 

To calculate the rhs of fTTTj) . it is very convenient to define the matrix- valued operator 
M to F by 



f dF 



(M(F)) 



dF 



dpjk dp jk 
dF 

dpkj 
dF 



{ d Pjj 

By (E2D, the rhs of (PU) is calculated to be 



(1 < j < k < m) 
(1 < k < j <m) 
(j = k = 1,2, - ■ ■ 



(22) 



m) 



d_ 
~dt 



F{l{t)) 



t=o 



E 

l<a<b<m 



dF 



dF 



dF 



dx 



ab 



(pMKb) + —(pMz'ab) \ + E q^(p) s : 



dy, 



ab 



r=l 



E {|£w5( s - + --) + 

l<a<»<TO k a0 



+ E O-(P) 
r=l u ^r 



E 

Ka<6<m 



<9F 



<9F , 4=J - 



9p a 



dp, 



«6 



= E {(^(^))^ + (^(F)UHL} + E(^(^))^; r 

l<a<b<m r=l 

= tr(.M(F)H'), (23) 

where the symbols, 9ft and 3 , stand for the real part and the imaginary part, respectively. 
To calculate the lhs of (TTTD. the introduction of the Hermitean matrix, 



= tf((gradF)(p))/i, (24) 
is of great use. Indeed, Eq. (1241) is put together with (JTU1) . (TT3"]) and (fl4l) to show 

m 7T a/' 



(((gradF)(p),H'))^ = 2 £ „ „ 

m m 



2 E «%{ E ^> 

j,k=l °3 "r" ^ a,6=l 



= 2 E ( E ^(/i t )ia)H / a6 = 2tr((/i^t )H / 

a, 6=1 j,k=l 

where = is the Hermitean matrix defined to be 



6j + 6i 



(j,k = 1,2, ...,m). 



(25) 



(26) 
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Since the Hermitean form, 

(S, S') G T p P m x T p P m ^ tr(StS') G C, (27) 
is well-known to be non-degenerate, we easily see, from ( 1231) and ( 1251) . that the equation 

M{F)=2hgtf + 2cl (28) 

has to hold true, where c in is the constant emerging from the trace-vanishing property 
of H' G T p P m (see (J7J)). The value of c will be determined below, soon. From (1241) and 
( |28]) . turns out to take the form 

djk = + k ) {(tfM(F)h) jk - 2c5 ifc } 

= - (etfM(F)h+ rfM(F)he-2cO) (j,k = 1,2, ■ ■■), (29) 
2 v /jfc 

where <5j& denotes the Kronecker delta. Equation ( |29l) is combined with (124]) and ( fTUl) 
to show 

(gradF)(p) = i(p>f(F) +A<(F)p) - cp. (30) 

We are now in a position to evaluate the constant c by taking the trace in both sides of 
f l30|) . By a simple calculation, the c is determined to be 

c = tr(pM(F)), (31) 
that leads us to 

(gradF)(p) = \(pM{F) + M(F)p) - (tr ( P M(F)))p. (32) 
To summarize, we have the following. 

Lemma 2.2 If a gradient system on the quantum information space (P m , ((•, -))® F ) is 
associated with a potential function F, it is governed by the equation of motion 

% = -UpM(F)+M(F)p) + (tr ( P M(F)))p, (33) 



dt 2 

where M.(F) is the matrix defined by |H|). 

Remark 2 In the case of the GS-NVNE for example, Lemma 2.2 is not so effective 
to derive its gradient vector since the negative von Neumann entropy is hardly written 
in terms of the entries of p. Indeed, in Uwano et al 2007, the gradient vector was not 
calculated directly on the QIS but was done through its lifting to STQM and projecting to 
the QIS; For such calculation, the geometric setting reviewed plays a central role. This 
could provide a good account for our review in subsection 2.1. 



3. The Karmarkar flow in the QIS 



This is the core section of this paper, where the Karmarkar flow is realized on the QIS. 
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3.1. The Karmarkar flow: Review 

Following Nakamura (1994), we review the Karmarkar flow for the canonical linear 
programming problem 

rn 

minimize ^ CjXj 

3=1 

m 

subject to Xj = 1, Xj > (j = 1,2,---, m) (34) 
i=i 

of unconstrained case (Karmarkar 1990), where CjS are given nonvanishing constants. 
We note here that no additional linear constraint specifying a subspace of the m — 1 
dimensional canonical simplex, 

m 

S = {x G R m | ^ = 1, x,- > (j = 1, 2, ■ ■ • , m)}, (35) 
i=i 

is taken into account in the unconstrained case. 

A continuous limit of the Karmarkar projective scaling algorithm gives rise to the 
system of differential equations 

The family of trajectories governed by (1361) is what we refer to as the Karmarkar flow in 
the present paper. The system of differential equations ( 1361) is brought into the following 
gradient-system form, according to Nakamura (1993): To be free from singularities, our 
discussion will be made on the regular part, 

m 

S = [x G R m | x i = !> x i > = !> 2 ' " " " > m )}> (37) 
i=i 

of the simplex S henceforth. With S, we endow the Riemannian metric 



- -Cjx) + xj(J2 c fc x fc) C? = !> 2, • • • , m). (36) 



((u,u>)) s x mp = Y<— (»,*■' ZTJ) (38) 

3=1 X l 

where the tangent space, denoted by T X S, of S at x G <S is given by 



T X S = {u e R m | u o = °}- ( 39 ) 

3=1 



On 5, the system of differential equations f[36|) admits the gradient-system form if we 
take the function, 

K (x) = ^x T Cx (x e S) (40) 
as the potential, where C is the diagonal matrix of the form 

C = diag(ci,c 2 , • ■ ■ ,c m ). (41) 
The gradient vector (grad k) (x) at x G S is defined to satisfy 

K (a(t)) Cu'eTj), (42) 



(((grad K )(x), M '))f mp = j f 



t=o 
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where a(t) with a<t<b (a <0<b) is a. curve in <S subject to 

= u' (43) 



t G [a, b] i-> a{t) G S, tr(0) = x, ^ 



t=o 



(cf. (|T71) with ( ITS!) ). By a straightforward calculation analogous to that for (grad _F)(p) 
in Sec. 2, we obtain 

m 

((grad k) (x)) j = Cjx] - x 5 ( E c k^t) (j = 1, 2, • • • , m), (44) 

fc=i 

which yields fl36|) . See Appendix B for a detail of the calculation. 

3.2. The gradient system on the QIS realizing the Karmarkar flow 

The gradient system on the QIS that we are seeking is constructed in what follows. 

3.2.1. The Riemannian structures of S and the QIS Let us consider the submanifold, 

m 

V={peP m | p = diag(0 1 ,...,O: E^ = !> e k >0 (fc=l,2,...,m)}, (45) 

fc=i 

of the QIS, which is easily seen to be diffeomorphic to S, the regular part of the canonical 
simplex S. Indeed, we can find the smooth one-to-one and onto map, 

/i : x G 5 i > diag (x u • • • , x m ) G V C P m . (46) 

Restricting the quantum SLD-Fisher metric ((■, -))® F of the QIS to the submanifold T>, 
we can make T> the Riemannian submanifold, whose metric will be denoted by ((•, -)) D 
henceforth. Namely, on expressing the tangent space of V at 6 (cf. ( TTOT) ) as the subspace, 

m 

T e V = {Ze M(m, m)\z = diag (G, • • • , Q, E Q = °}> ( 47 ) 

of T@P m , ((•, -)) D i s defined to satisfy 

((Z, Z'))g = ((Z, Z'))g F (Z, Z' G T e P C T e P m ). (48) 
We show the following. 

Lemma 3.1 The map \i, defined by fcjfty , of S to T> is isometric; the identity, 

(CM**), )M«0 ))?(,) = «'))f mp («, «' e (49) 

ZioZds irue, where i/ie differential of the map fx at x defined by 

f**A u ') = ~J1 M^OO) = dia g K> • • • > M m) ( M ' G T x 5) (50) 
at t=0 

with (SW and (g^). 

Proo/: Equation (TJHD is put together with Eqs. ([ID])- (HID and flSDD to yield 



This completes the proof. 
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3.2.2. Construction of the gradient system From Lemma EH1 we learn the coincidence 
of the Riemannian structures of the regular part, S, of the canonical simplex for the 
Karmarkar flow, and of the submanifold, T>, of the QIS. Accordingly, in order to find a 
gradient system realizing the Karmarkar flow on the QIS, we naturally come to seek a 
function K(p) on the QIS whose restriction to T> coincides with the potential k(x) for 
the Karmarkar flow through the map p: We take K(p) to be 

K(p) = 1 -tr(Cp 2 ), (52) 



2 

where C is the diagonal matrix given in (1411) . The entries of the matrix M.(K) given 
by (|22p with K in place of F are calculated to be 



(M(K)) jk = = hcj + c k )p kj = hcj + c k )p jk (j, k = 1, • • • , m), (53) 
which is brought into the form 

d -£ = \(C P+( C). (54) 
Equation (l54l) is put together with (132]) to show 



(gradK)(p) = ](p 2 C + 2pCp + Cp 2 ) - (tr ( P Cp))p, (55) 



4 

so that we have the following. 

Lemma 3.2 The gradient system on the QIS associated with the potential K(p) is 
governed by the equation of motion, 

f = ~\(P 2 C + 2pCp + Cp 2 ) + (tr (pCp))p. (56) 

We are at the final stage to see how the gradient system on the QIS with K realize the 
Karmarkar flow. Recalling Eq. (1551) . we immediately obtain 

(gradK)(6) = CQ 2 - (tr (CQ 2 ))e G T @ V C T e P m (57) 

(see ( 1471) for TqV), which enables us to restrict the equation of motion (J56"j) to the 
submanifold T> of the QIS isometric to S. The restriction indeed gives rise to the 
system of differential equations, 

d9i 



'3 



dt 



-c 



•A 2 + M£ c ^) (J = l,2,---,m), (5* 



k=l 

on T>, which is evidently identical with the Karmarkar flow ( 1361) with 9 in place of x. In 
conclusion, we have the following. 

Theorem 3.3 The gradient system on the QIS associated with the potential K(p) (GS- 
QIS-K) realizes the Karmarkar flow on the submanifold T> of the QIS. 
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4. Concluding remarks 

We have successfully constructed the gradient system (GS-QIS-K) on the QIS which 
realizes the Karmarkar flow on the submanifold T>. A key to the success is the isometry 
of the underlying Riemannian manifold S for the Karmarkar flow and the Riemannian 
submanifold T> of the QIS, that is presented in Lemma [3.1[ 

Through the study leading us to Lemma 13.11 on the Riemannian structures of S 
and T>, a clear account for the encounter with the GS-NVNE as a generalization of 
the GS-MD is obtained: The Riemannian metric and the potential for the GS-MD are, 
respectively, equal to those for the GS-NVNE restricted on D up to a common multiplier. 

Integrability of the GS-QIS-K is an open question. We wish to recall that, in 
the case of the GS-NVNE (Uwano 2006), the invariance of the negative von Neumann 
entropy chosen as the potential under the U(m) action, p i— > hph* (h G U(m)), works 
effectively to show the integraility in the sense that the GS-NVNE allows a sufficient 
number of mutually independent integrals of motion. The U(m) invariance of the 
potential K(p) does not hold true, however, so that we have little expectation of the 
integrability for the GS-QIS-K. If we find a U(m)-invariant potential whose restriction to 
T> realize k(x), the gradient system with that potential on the QIS would be integrable 
and, further, would admit a double-Lax bracket representation like in the case of the 
GS-NVNE. 

In the case that the matrix C in the potential K(p) is taken to be C = 21, 
the GS-QIS-K has two other particular features: One is that it turns out to realize 
not only the Karmarkar flow but also the flow solving the eigenvalue problem of anti- 
Hermitean matrices in view of Nakamura (1992, 1993). The other is that the potential 
K(p) = tr(p 2 ) is called the purity whose logarithm provides the minus of the quantum 
Renyi potential log(trp 9 )/(l — q) with q = 2; the larger the purity of a quantum state 
is, the larger its Hilbert-Schmidt distance from the maximally mixed state. (Bengtsson 
and Zyczkowski 2006). Due to the second feature, the potential K(p) with C = 21 could 
be interpreted to be an object in quantum physics. A paper dealing with the case of 
C = 21 is in preparation. 

The averaged learning equation of Hebb-type dealt with in Nakamura (1994) is 
a current target of the authors: Through another geometric trick, we have recently 
succeeded to find its generalization on the QIS, which will be reported in other paper 
(in preparation). 
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Appendix A. The Riemannian metric ((•, -)) R 
Appendix A.l. Geometry of ir^(P m ) 

We start with introducing the natural Riemannian metric of Mi(2 n , m) (d 7r~ 1 (P m )). 
On regarding M(2",m) = C 2 " as the 2 n+1 -dimensional Euclidean space, the real part of 
the Hermitean inner product of M(2 n , m) given by ([1]) can provide the Euclidean metric 

(X,X')f = -^tr(XtX' + (XtXO) 
2m 

(X, X' G M(2 n , m) T$M(2 n , m), $ G M(2 n , m)), (A.l) 

where T^C 2 ™ denotes the tangent space of M(2 n ,m) at $ G M(2 n ,m). The tangent 
space, T$M!(2 n ,m), of Mi(2 n ,m) at $ G Mi(2",m) is thereby defined to be 

T$Mi(2 n , m) = {X G M(2 n , m) | tr(X f $ + $ f X) = 0}, (A.2) 

on looking M!(2 n ,m) upon as a submanifold of M(2 n ,m). Through the restriction of 
T$M(2",m) to T M 1 (2 n ,m) ($ G Mi(2 n ,m)), the Euclidean metric (•, -) E of M(2 n ,m) 
is brought to the Riemannian metric of of M!(2 n ,m), 

(X, X')f = — tr {pX' + (XI X>)) (X,X' GT $ M 1 (2 n ,m), $GM 1 (2 n ,m)). (A.3) 
2m 

On account that T^n m \P m ) = T^M 1 (2 n ,m) for $ G ^(PJ (c M 1 (2 n ,m)), the (•, -) R 
becomes the Rimennian metric of the inverse image n^(P m ) of P m by Tr m if restricted. 

Appendix A.2. The horizontal lift 

We introduce the horizontal lift of tangent vectors of T p P m to T^,n^-(P m ) (7r m ($) = 
p G P m ) as follows: According to the U(2") action ([5]), let us consider the orthogonal 
direct-sum decomposition 

T $ M 1 (2 n ,m) = Ver($)© ± Hor($) ($ G 7r~ 1 (P m )) (A.4) 

for the metric (•,•)$ with 

Ver($) = {X G T^\P m ) \X = r ] <S>, V e u(T)} ($ G tt^PJ) (A.5) 
and 

Hor($) = {X G T^-\P m ) | $Xt - X$t = 2 ,, 2 n} ($ G tt" 1 ^)), (A.6) 

where u(2 n ) denotes the set of all the anti-Hermitean matrices of degree-2 n and C>2«.2™ 
does the null matrix of degree-2 n (see also Prop. 4 in Uwano et al 2007). The horizontal 
lift, denoted by ^$(S), of H G T p P m is then defined to be the unique tangent vector at 
$ G vim. 1 (p) subject to 

ir m ^(U(Z)) = E and £$(H) G Hor($). (A.7) 

The 7r m!|<j( j ) is the differential of the map ir m at $ G P m , which is defined to be 

d 



7r m*,$(A A ) — — 



7T m ( C (t)) (XGT $ P m ), (Ai 
t=0 
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where c(t) with a < t < b (a < < b) is a curve in Mi(2 n , m) satisfying 

dc 



t G [a, b] i — > c(t) G Mi(2 n , m), c(0) = $, 



= X. (A.9) 

t=o 



Under the expressions ffTUl) and (fTTj) of p G P TO and H G T p P m , the singular- value 
decomposition (see Rao and Mitra 1971), 

^sf^)^ ( 3 <?GU(2")) with v / e = diag( v /^,---, v / t)' ( A - 10 ) 

\ < ^2 n -m,m ) 

of $ G vr~ 1 (p) works effectively to obtain the horizontal lift. Indeed, on putting flA.101) 
with (TlOl) and ( TTTT) together, the horizontal lift of H G T p P m is given by 

4(2) = ( + ^ x)) ) (A.ll) 

where ct^Q>(x) stands for the m x m anti-Hermitean matrices uniquely determined by 
e- 1 a V e( X ) + a^We- 1 = -O" 1 * + X®' 1 (A.12) 



(see Uwano e£ a/ 2007). By a straightforward calculation, we have 

' 0j - e_ 

9i + 9, 



( a Ve(x)) ik =(^—Y)x jk (j,k= 1,2,-.., m). (A.13) 



Appendix A. 3. Defining ((•, 



On using the Riemannian metric (•, -) R of Mi(2 n ,m) and the horizontal lift £<&(■), the 
Riemannian metric ((•, .)) R of P m is defined to be 



(U(E), (P e ^m, 2, 5 G T p P m ) (A.14) 

where $ G tt" 1 ^) can be chose arbitrarily. Equations (1A.3I) . (jA.llD . (1A.13I) and (1A.14I) 
are put together to yield ffT5]) . 



Appendix B. The gradient vectors 



We calculate the gradient vector (grad/t)(x) at x G <S of k according to (f42|) . On taking 
(T4"3"j) for the curve a{t) into account, the rhs of (T4"2"j) is calculated to be 

Further, Eq. ( 1381) is combined with the lhs of ( f42i) to show 



t=o i=i ut j=i 



rn 



(((grad«)(x),« / ))* mp = E — • ( B - 2 ) 



j=i x i 



Equations (IB.lj) and (IB. 21) therefore yields 
((grad/c)(a:)) i 



CjXj + c, (B.3) 
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where c is a constant common in j (= 1,2, •• -,m) emerging from the constraint 
YIj=i u 'j — to u' G T X S (cf. the constant c in Eq. fl28l) ). Indeed, to fulfil the condition 
J2]Li ((grad k)(x))j = for (grad/t)(x) to be in T X S, the constant c turns out to satisfy 

m m m m 

= £ ((S rad K ) ( x ))j = £ c i x ? + 5 £ x i = £ c i x ? + 5 ' ( B - 4 ) 
i=i i=i i=i i=i 

which is solved to be 

m 

i=i 

In the sequel, Eqs. (1B.3j) is put together with (1B.5[) to show 
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